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1 Introduction 



Application of the interacting boson model (IBM) [|1| to deformed nuclei has been 
fraught with conceptual and technical difficulties from the early days. Truncation of 
the model space to s and d bosons created a lot of controversy and was challenged 
on both microscopic ^ and phenomenological |^ grounds. When the dust finally 
settled, it was generally accepted that g bosons were both necessary and sufficient 
for microscopic foundation |Q as well as phenomenological description of deformed 
nuclei in the IBM ||^ (see also P, §| for recent reviews). The sdg-IBM, however, 
does not share many of the simplifying features of the sd-IBM which made it so 
popular. To start with, the symmetry limits of the dynamical group SU(15) are not 
associated with actual spectra so there is little help from group theory. Secondly, 
due to the large basis spaces, exact diagonalization of sdg-IBM Hamiltonians is not 
possible for deformed nuclei. As a result, progress in sdg-IBM calculations has been 
rather slow, and due to various approximations involved, a satisfactory description 
of both low-lying band structures and high-spin states in deformed nuclei is still 
missing. 

Numerical calculations in the sdg-IBM were initially performed by coupling a 
single g boson to an sd boson core . In another approach, a Hamiltonian consist- 
ing of various SU(3) tensor operators was diagonalized in a truncated SU(3) basis 
1^. At present, full basis sdg-IBM calculations are possible only for vibrational and 
transitional nuclei with boson numbers < 10 [0. The accuracy of the truncated 
space calculations for deformed nuclei with A^ > 12 has not been explored so far. 
One of the purposes of this paper is to examine the convergence properties of various 
matrix elements (m.e.) with the number of g bosons allowed, and to establish some 
criteria for the accuracy of truncated sdg-IBM calculations. It should be clear that 
due to truncation, the validity of numerical diagonalization results in the sdg-IBM 
is limited to low-lying states. Besides, such calculations are time consuming and 
therefore are not well suited to explore the effect of the multitude of parameters in 



the sdg-IBM. The 1/A^ expansion method |Tl|, which is based on angular momentum 



projected mean field theory, could alleviate both problems. The available analytic 
formulas for various physical quantities allow systematic study of model parame- 



ters, and fast and efficient analysis of data [12|. In addition, the 1/N expansion 
has recently been extended to higher orders using computer algebra [0, thus en- 
abling accurate descriptions of high-spin states within the sdg-IBM. Currently the 
expansion formalism covers ground- and single-phonon bands. Calculations 
for multi-phonon bands are rather complicated and have not been performed yet. 
Therefore, at this stage, one has to rely on both analytical and numerical methods 
to obtain a complete description of deformed nuclei 

In this paper, we ffist discuss the recent developments in the numerical diag- 
onalization and the 1/N expansion methods. We then present a systematic study 
of deformed nuclei in the sdg-IBM employing a Hamiltonian consisting of one-body 
energies and quadrupole and hexadecapole interactions. Application of the results 
is focused on the high-spin states in rare-earth and actinide nuclei as their descrip- 
tion in the IBM has been a source of criticism [^] which has not been adequately 
addressed so far. In order to constrain the model parameters properly, both the 
high-spin data and the low-lying band structures are described simultaneously. 

2 Choice of Hamiltonian 

The consistent-Q formalism (CQF) |^ has been the standard choice for description 
of deformed nuclei in the sd-IBM. The CQF Hamiltonian consists of the dipole and 
quadrupole interactions 

H = K^L-L-K2Q-Q, (1) 
where the quadrupole operator is given by 

Here brackets denote tensor coupling of the boson operators and 6^^ = (— 
For consistency, the same quadrupole operator is used in the E2 transition operator, 
T(E2)= 62(5, as in the Hamiltonian. The CQF has been especially successful in 
reproducing the low- lying band structures and the E2 transitions among them 0. 
A more complete description of deformed nuclei, including high-spin states and 
hexadecapole bands, however, requires extension of the model space to sdg bosons. 
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A minimal extension of the CQF Hamiltonian to the sdg-IBM can be achieved 
by including the g-boson energy term, Egfig in the Hamiltonian (|l|), and modifying 
the quadrupole operator to 

Q = [,tj+ rftg](2) + q^^ld^df) + q^^ld^~g + g^f) + q,,y~gf). (3) 

We shall refer to this minimal extension as the CQF below. A study of high-spin 
states in the sdg-IBM using the CQF Hamiltonian has indicated that the energy 
surface remains too rigid, and inclusion of the d-boson energy term, Edfid, is essen- 
tial in order to reproduce the spin dependence of moment of inertia (MOI) and E2 



transitions in the ground band |]T3[. The success of this pairing plus quadrupole type 
of Hamiltonian, however, does not extend to the side bands which are more sensi- 
tive to interference from the hexadecapole interaction. Thus, for a comprehensive 
description of deformed nuclei one needs to employ the Hamiltonian 

H = ejid + eghg - K2Q ■ Q - ■ T^, (4) 

where the hexadecapole operator is given by 

T4 = [s^g + g^sp + h,2[d^dp + h,,[d^g + ^?tj](4) + h,,[g^gp. (5) 

Note that we have deliberately left out the dipole interaction in eq. @ as it is at the 
root of the rigid MOI problem in the IBM. This Hamiltonian contains 10 parameters, 
namely, the one-body energies Sd and e^, the multipole interaction strengths K2 and 
K4, and the quadrupole and hexadecapole parameters qji and hji. In a systematic 
study covering many nuclei, it is desirable to have a smaller set of free parameters. 
To achieve this goal, we adapt a similar strategy as in a previous study of deformed 
nuclei [0. The quadrupole parameters {522, ?24, 944} are scaled from their SU(3) 
values with a single factor q as suggested by microscopies [0 , while the hexadecapole 
parameters hji are determined from those of qji through the commutation condition, 
\h, q\ = 0, which ensures that the quadrupole and the hexadecapole mean fields are 
coherent 

h22 = q24, h24. = qu, ^^44 = ^24 + - ^22^44 - l)/g24 • (6) 

Here qji = {j0l0\20)qji and hji = (jO/0|40)/ij;. This reduction of parameters from 
10 to 5 is obtained at the expense of detailed description of quadrupole and hex- 



adecapole operators. Since information, especially on the latter, is rather patchy, 
this will not cause any problems except in a few isolated cases. In the calcula- 
tion of E2 and E4 transitions, we shall use the consistent operators, T(E2)= e2Q, 
r(E4)= 64X4, so that, apart from the effective charges 62 and 64, no new parameters 
are introduced. 

In concluding this section, we present an alternative parametrization of the 
Hamiltonian (Q) which is more convenient in implementation of the expan- 
sion formulas. This involves factoring out the energy scale and the leading order 

dependence from the energy expressions. Since the quadrupole interaction is 
dominant, a suitable choice for such a set of dimensionless parameters is given by 

r]i = ei/NK2, Ck = /tfc//^2, (7) 
where Z = 0, 2, 4, . . . correspond to the subscripts s,d,g, . . .. 

3 Numerical diagonalization 

A computer code which can diagonalize arbitrary sdg-IBM Hamiltonians in full space 



has been available for some time [|T^], but due to excessive memory requirements it 



had only limited applications to transitional nuclei |jTO[ . Recently this code has been 



modified to run on supercomputers improving its applicability [0. Nevertheless, 
exact diagonalization for deformed nuclei with > 10 remains elusive due to the 
large basis space, and truncation of the model space is still necessary. Here we use 
the modified code to study the accuracy of the truncated space calculations in the 
sdg-IBM. Numerical results will also be used in demonstrating the accuracy of the 
1/A^ expansion formulas in sect. ^, and in application to actual spectra in sect. ^. 

Since g bosons are relatively weakly coupled, a natural parameter in truncating 
the basis space is the maximum number of g bosons allowed, rigmax- In figs. 1-2, we 
examine the convergence properties of some key observables as a function of rigj^^x- 
For this purpose we use a typical deformed Hamiltonian with parameters (as defined 
in eq. (0)), H2 = -20 keV, ^4 = 0, g = 0.5, r]d = 1.5, % = 4.5, and A^ = 10. Fig. 1 
shows the effect of the truncation on low-lying band structure: a) band excitation 
energies, b) E2 transitions, and c) E4 transitions. The rigmax = 1 calculations are 



off by about 10-20% (mostly overestimated but underestimated in a few cases), and 
hence they are not very rehable. As expected, the hexadecapole bands take longer 
to converge compared to the (3 and 7 bands, the worst case being the (3' band. 
Nevertheless, convergence to accuracy of a few percent is obtained in almost all 
cases for rigmax = 3. In fig. 2, we present a similar study for the high-spin states in 
the ground band: a) excitation energies, and b) E2 transitions. At spins L ~ 2A^, 
the rig max = 1 calculations are off by about 20-30% which will get even worse with 
increasing spin. The rigmax = 3 results, on the other hand, provide a reasonably 
accurate picture up to spins L ~ 2A^. Beyond that, g bosons start dominating the 
wave functions, and any truncation is likely to lead to substantial errors. 

The above results suggest that diagonalization of the sdg-IBM Hamiltonians in 
a model space truncated to rigmax = N/3 bosons will give a reliable description of 
states with spins L < 2N. This extends the applicability of the SDGBOSON code 
to = 14 which covers roughly half of the deformed nuclei. It should nevertheless 
be emphasized that these computations are expensive, time consuming, and cer- 
tainly not the best way to deal with the sdg-IBM problems. In the next section, 
we introduce the 1/A^ expansion which circumvents the shortcomings of numerical 
diagonalization. 

4 expansion formalism for high-spin states 



The 1/A^ expansion formalism |Tl| was developed as a response to difficulties in 
performing calculations in the sdg-IBM due to the inadequacy of group theoretical 
techniques and the large basis space problem in numerical diagonalization. It is 
based on angular momentum projected mean field theory and leads to analytic 
expressions for various physical quantities of interest. Initially, the 1/A^ calculations 
were carried out to order 1/A^^ which is quite sufficient for low-lying states. An 
accurate description of high-spin states, on the other hand, requires inclusion of 
terms up to order 1/A^^ which is not suitable for hand calculation. This difficulty 
has finally been overcome through the use of computer algebra. A brief description 



of the method for the ground band was given previously [T^. Here we provide 
details of the extended calculations including the results for the excited bands and 
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the electromagnetic transitions. 



4.1 Ground band 

We consider a general formulation of the IBM as this allows an elegant derivation 
of the formulas by fully exploiting the angular momentum algebra. Thus we 
introduce the boson creation and annihilation operators b]^, bi^ with / = 0, 2, 4, . . .. 
where bo = s, 62 = d, 64 = g, etc. In order to keep the variational problem 
to a manageable size, it is necessary to assume that the boson system is axially 
symmetric, and hence K is a good quantum number. From comparison with the 
exact diagonalization results, this assumption will be seen to hold to a very good 
degree. The ground band then can be written as a condensate of intrinsic bosons as 

\<P,) = {N\)-y'{bY\0), b^ = T.^ibl (8) 

I 

where xi are the normalized boson mean fields, i.e. x.x = 1, x = {xq, X2, X4, . . .). In 
the classical limit of the IBM, the mean fields are associated with the deformation 
parameters of the system [113. For a given Hamiltonian H, they are determined 



from {H)l by variation after projection (VAP). 

The Hamiltonians in eqs. can be written in the generalized form as 

2/ max 

H = E^ini - E '^'^T^'^ ■ ^^'^ = j:blk,, T^'^ = j:tUb]bif\ (9) 

I k=0 fJ, jl 

where the parameters have the obvious correspondence, £2 = ^d, ^4 = £g, hji = Qji, 
t^ji = hji. This general form has the advantage that, to evaluate the expectation 
value of H, one needs to perform the calculation for a generic number operator ni 
and a multipole interaction T^'^^ ■ T^'^\ The expectation value of a scalar operator O 
in the ground band (|^), with angular momentum projection, is given by 

or _|_ 1 r 

{0)l = 2NW{<p^, L) I ^/?sin/34(/5)(0|6^Oe-^^^^(6t)^|0). (10) 



Here, the normalization, M{(j)g, L), follows from eq. (10) upon substituting the iden- 



tity operator for O. Algebraic manipulations in eq. ([10| ) are most easily carried out 



using boson calculus and angular momentum algebra techniques (see ref. [11] for a 
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pedagogical treatment). For the number operator, one obtains 

{ni)L = E(L0/00|/0)2F(iV - 1, /), (11) 

where F{N, L) denotes the reduced normalization integral 

F{N,L)=Af{(Pg,L)/{2L + l). (12) 
Eq. ([IlD is exact and highlights the essential role played by the normalization inte- 



gral. In the original papers a Gaussian approximation was used in the evalu- 
ation of F{N,L) which limited the accuracy of m.e. to order 1/N'^. This difficulty 



has been overcome recently using the computer algebra software Mathematica |T9| 



By exploiting the symmetries of the boson system, the normalization integral is cast 



into a system of linear equations which is solved with the help of Mathematica pO |. 
The result is a double expansion in 1/A^ and L = L{L + 1) given by 

F(N, L) = 4? y , y anmL"^. (13) 

The coefficients a„m in eq. (|T3|) are given in terms of polynomials of the moments 
of xf 

a^ = Y^r+^xl (14) 
I 

and a is defined as a = oq. A list of a„m up to sixth order is given in ref. [^. The 
knowledge of F{N,L), in principle, allows evaluation of m.e. to arbitrary orders in 
1/A^. As will be seen in the applications, a correct description of MOI at high-spins 
requires inclusion of terms of order L^/N^. Evaluation of eq. (|ll]) to such high orders 
is too difficult to perform by hand but becomes manageable using computer algebra. 

Before presenting the final results, it will be useful to comment on the general 
form of the m.e. of a /c-body operator O, and illustrate the concept of layers in the 
1/A^ expansion 



+(^no3o + ••■)+■■•}. (15) 

The expansion coefficients Onm in eq. (|T3|) involve various quadratic forms of the 
mean fields xi corresponding to the single-boson m.e. of O and its moments. The 
explicit form is given to facilitate the illustration of layers. Notice that the i co- 
efficients Onm in the i'th column have n + m = i — 1 constant, and are referred 
as the layer "z — 1". The leading term in eq. (15) thus forms the zeroth layer. 



(This name is appropriate since calculations in the intrinsic frame give the same 
result independent of projection.) In simple terms, the layer in an expansion is 
given by the maximum power of L. There is a close connection between the lay- 
ers in the m.e. ([T5|) and the normalization coefficients anm in eq. (|13D, namely, in 
order to calculate the m.e. up to the z'th layer, one needs to know the coefficients 
{arm, Oinn-1, ■ ■ ■ , «nn-i+i, n = 1, 2i} . This is Very useful in higher order calculations 
as it restricts the number of terms in the expansion, cutting down the amount of al- 
gebra. To make this point clear, we note that eq. ([I5D shows all the terms in the third 
layer whereas a complete calculation to order would require 6 more terms be- 
longing to the fourth, fifth and sixth layers. As can be seen from eq. (|T^) below, the 
complexity of the coefficients Onm increases "exponentially" with layers, and each 
of the extra terms would lead to expressions pages long. From a practical point 
of view, such accuracy is never required. The only term of any consequence 

is L^/N^ which is included in the third layer. The rest are completely negligible. 
Hence use of layers is a more sensible approach than a complete calculation to a 
given order in 

With these considerations, we present the result of the Mathematica evaluation 
of the one-body m.e. ( pT]) to the third layer 

{hi)L = Nx^{l + -L(^a - Z) + + ai/2 + (1 - a,/a)l + F/2) 

(a + 2a^ — 7ai/3 — aai + ha\/Aa — 02/8 



{aNf 

+ (-1 - 2a + 2ai + 7ai/2a - ^al/2a^ + a2/2a)l 



+ (-7/6 - a + 5ai/4a)P - P/6) 



L 
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-a 



2>{aN) 



I) + ^(2a + - 2ai + (-2 - 2a + 3ai/a)/ - F 

— ^f-3a - 12a2 - 4a^ + 21ai/2 + llaai - 15a?/2a + 802/2 
[aNy ^ 

+ (3 + 12a + 4a^ - 33ai/2 - 14ai/a + 25ai/2a^ - 2a2/a)r 
+ (7/2 + lla/2 - hai/a)P + f /2) 

' -a - 2a2 + 3ai/2 + (1 + 2a - 2ai/a)r+ P/2) 

f^(4a + 21a^ + 14a^ - 16ai - 51aai/2 + 13a^/a - 2a2 

+ (-4 - 21a - 14a2 + 34ai + 20ai/a - ?,9al/2a^ + 5a2/2a)r 
+ (-4 - 17a/2 + 13ai/2a)P - P/2^ 

a - 6a^ - 6a^ + 25ai/6 + 9aai - 15ai/4a + 5a2/12 
+ (1 + 6a + Qo? - 45ai/4 - 5ai/a + 21ai/4a^ - a2/2a)l 
+ (5/6 + 9a/4 - 3ai/2a)/^ + P/I2] }, (16) 



2{aNy 
"aN 



where a^ is defined in eq. ([T^) . Eq. ([T6|) can be checked against two results: i) 
it satisfies the number conservation, i.e. J2i{fii)L = N, and ii) it reproduces the 
analytic formulas available in the SU(3) limit [Q]. 

A similar calculation for the multipole interaction yields the intermediate result 

(T« . T«), = ^^^|^{E %ff ULOmiOmN - 1, /) 



+ iN-l) J2 tkjitkfi'XjXix,,xv{jOj'0\JO){lOl'0\JO) 
jlj'l'J 



I' I k 



J2{L0J0\I0)^F{N -2,1) \. 



(17) 



Again this is exact and can be evaluated to any order using Mathematica. The third 
layer result is given by 



(aA^) 



(aA^) 



^((-2a + ai)f/fc + (1 - a - aJa)Uki + Uk2/2 + a^Ck - aCki) 

^f(2a + 2a2 - 14ai/3 - aai + 5a2/2a - 2a2/3)f/fc 
/V w V 

-1 + a - ai/2 + 7ai/2a - 5a^/2a^ + a2/2a)Uki 
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{aNy 



+ (-7/6 + 5ai/4a)[/fc2-f/fc3/6 
+ {-a^ + aai/2)Ck + {a - ai)Cki + aCk2/2] 

2aUk + Uki 

^((4a + 2a^ - 4ai)Uk + (-2 + a + 3ai/a)Uki - Uk2 - a^Ck + aCki) 



aN 



(aN) 



— r((— 6a — 16a^ — 4a'^ + 21ai + 15aai — \ha\/a + 802)?/^ 



2{aNy- 



+ (3 + 2a - 2a^ - 4ai - 14ai/a + 25ai/2a^ - 2a2/a)Uki 
+ (7/2 + 2a - 5ai/a)f/fc2 + Uk^/2 

+ {2a^ + 2a^ - 2aai)CA: + (-2a - 2a^ + 'iai)Cki - aCk2 

{-2a - 2a^ + 3ai)Uk + (1 - 2ai/a)Uki + Uk2/2 
1 



aN 



(8a + 30a^ + 14a^ - 32ai - 37aai + 26ai/a - Aa2)Uk 



?>{aNy 



+ (-4 - 8a + 2a^ + 29ai/2 + 20ai/a - ?,9al/2a^ + 5a2/2a)f/fci 
+ (-4 - 9a/2 + 13ai/2a)f/fc2 - t/fcs/S 
+ (-a2 - 2a^ + 3aai/2)Ck + (a + 2a^ - 2ai)Cfci + aC;.,2/2; 

-8a - 36a2 - 24a^ + lOOai/3 + 54aai - 30a^/a + 10a2/3)[/fc 



+ (4 + 12a - 24ai - 20ai/a + 21ai/a^ - 2a2/a)U, 



kl 



+ (10/3 + 6a - Qai/a)Uk2 + f/fcs/s] }• 



(18) 



Here the quadratic forms Ckn arise from normal ordering and simulate an effective 
one-body term 

Ckn = {2k + l)Y,r{tk,iXif/{2l + 1), (19) 
while Ukn represent the genuine two-boson interaction 

Ukn= E 7"(j0j'0|/0)(/0r0|/0)| ^'^ ^ "^tkjitkfi'XjXiXfXv. (20) 

For a given multipole, these sums can be evaluated in closed form using Mathemat- 
ica. For the quadrupole and hexadecapole interactions, the first four terms needed 
in eq. (|18D are given by 

U2 = A\ 
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U2i = (2A - 3A)A, 

U22 = {2A2 - 24Ai + \%A)A + (An - + 7Ai)Ai + (An - ^2)712, 

\J2z = (2A3 - 36^2 - 18An + 240^1 - 144A)A 

+ (3^21 - 3A3 + 56^2 + IQAn - 194Ai)Ai/2 
+ (llA?i + 14AnA2 - 25A2)/12 + {A3 - A2i){A2 - An)/4, 
= B\ 

Uii = {2Bi - 10B)B/2, 

Ui2 = 4:Bl + (2fin - 40fii + 205)5 + {B2 - Bn - 205i + 1805)Vl80 

^743 = (2^3 - 120^2 - 60Bn + 24805i - 44005)5 

+(3521 - 353 + 22452 + 585n - 27565i/3)5i/9 

+ (85i\ + 145n52 - 1152)/45 + (53 - 52i)(52 - 5n)/60 (21) 

Here the quadratic forms A^n and Bmn in (|2T|) are defined as 

Amn = Ej"^""(j0^0|20)t2,7X,x;, 5„„ = Y.j^'r {jOlO\AO)UjiXjXi, (22) 
jl jl 

and correspond to various moments of the single-boson m.e. of the quadrupole and 
hexadecapole operators, (note that the zero subscripts are suppressed for conve- 
nience). The quadrupole m.e. given by ( [I8| - |2TD reproduces the well known Casimir 
eigenvalues in the SU(3) limit, hence also passes the SU(3) test. 

The analytic expressions presented above are already rather long. If for any 
reason, the next layer results should be required, the expressions would grow to 
pages long, and the analytical 1/A^ calculations might not be very practical. In such 



cases, numerical evaluation of the m.e. (pA|jT7| ), as described in Appendix A, may 



be preferable. Although this would increase the computation time appreciably, it 
has the advantage that the calculations are done exactly to all orders in 1/A^. 

4.2 Variation after projection 

The energy expression derived in the last subsection is rather lengthy, and in dis- 
cussing the variational problem, it will be more convenient to express it in a compact 
form. Thus, using the parametrization in eq. (|^), we rewrite the ground band energy 
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as 



(23) 



where the coefficients Enm can be read off from eqs. ([T6| - p2D . For example the 
leading order is given by 

.2 . ^ .2 X 5 .2 



00 



x.x 



x.x 



x.x 



(24) 



with A and B defined in eq. (^2]), and we have restored the normalization factors 
x.x as a precursor to variation. The minimum of the ground energy is obtained from 



dE^L/dxi = 0, 1 = 0,2,4, 
which can be solved algebraically using the ansatz 

E^nm ( ^ \"' 
Mm \ Jj2) ■ 



(25) 



(26) 



The use of layers again simplifies solution of the variational equations. For the 
leading order (zeroth layer), one has the usual Hartree-Bose equations 

OEqq 



dxi 



0, 



(27) 



XOO 



which are a system of coupled non-linear equations, and they are solved numerically 
by iteration Having determined xqo, the first layer mean fields xqi and xio are 
then obtained by solving the respective sets of equations 



dE, 



00 



dxi 

OEqq 



1 dE, 



01 



xoo+xoi/Af 



dxi 



N dxi 
L dE^o 



XOO 



m dxi 



(2J 



XOO 



xoo+xioL/Af2 

Upon substituting the mean fields in derivatives in (^Sf), the leading order vanishes 



by virtue of the Hartree-Bose eqs. (27), and the next order leads to sets of linear 
equations for xqi and xio, that can be easily solved using Mathematica. The Hartree- 
Bose condition also ensures that when the first layer mean fields are substituted in 
the energy expression, the correction to the first layer exactly vanishes. So they only 
contribute to the second and higher layers [O]. This holds in general for all layers. 
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Thus for the third layer expansion considered here, one needs at most the second 
layer mean fields xo2, xn and X20 which are obtained from 



00 



dxi 



1 dE, 



01 



xoo+xoi/Af+xo2/iV2 



1 dE, 



02 



dE, 



00 



dxi 



N dxi 
L dEio 



xoo+xoi/A^ 



dxi 



xoo 



dE, 



00 



dxi 



xoo+xioL/Af2+X2oi2/A''' 
xoo +X01 /Af +X10 L/Ar2 L/Ar3 



dxi 

1 9^01 



L2 dE~ 



20 



xoo+xioL/A^^ 



xoo 



dxi 
L dEn 



xoo+xioi/Af^ 



dx, 



iV3 



xoo+xoi/A'' 

(29) 



xoo 



Again these sets of linear equations can be solved using Mathematica. We refrain 
from presenting these rather bulky results for the first and second layer mean fields 
here because, in the absence of analytical solutions for the zeroth layer, they are not 
very illuminating. Upon substituting eq. (^) in (^), one obtains the variational 
corrections introduced by the higher order mean fields in the ground band energies. 
These lengthy analytic expressions contribute only to the second and higher layers 
and will not be shown here. All these results, together with other 1/N expansion 
formulas, are nevertheless available in the form of a Fortran code [^. Finally, if 
one is interested only in practical applications of the results to high-spin states. 



one can determine the minimum directly from the energy expression ([23| ) using the 
numerical simplex method, and thereby avoid the complexities introduced by the 
higher order terms in the solution of the variational problem. 

4.3 Single-phonon bands 

Most of the high-spin data, as well as their theoretical analysis, are concentrated on 
the yrast bands (ground or two-quasiparticle) . While relying solely on the yrast data 
may be tolerated for microscopic models, it could easily lead to misleading results in 
phenomenological models. For this reason, inclusion of single-phonon bands in the 
analysis of high-spin data is highly desirable in phenomenological approaches. As 
will be seen in the applications, there are substantial high-spin data for the 7 bands, 
which can be singled out among the single-phonon bands in this respect. Therefore, 
we consider here the 7 band as an example of a single-phonon band calculation. 
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Energy expressions for the other bands can be derived in a similar fashion. 

The single-phonon bands are obtained from the ground band by acting with 
the other intrinsic boson operators b}^ = J2i Ximb\m^ and then orthogonahzing the 
resulting bands. For example, the 7 band intrinsic state is given by 

10,) = iV - 1) + -^^=^e,(6l)^|0„ iV - 2), (30) 

In this trial state, the mean fields for the ground band are already established in 
the last subsection, and those for hi are determined from the spurious K = 1 band 
as [1TI| , 

xii = [l/aY/^xi. (31) 

Thus, only the 7 band mean fields, X12, are to be determined by VAP. The coefficient 
in eq. (|30|) follows from the orthogonality condition {L^\Lg) = as 



5^(L2/-2|/0)(L0/0|/0) 
II 



xiXi2F{N -1,1) 

^ (32) 



+e,5]a;,x,ix,vx,vi(jO/0|/0)(jl/l|/2)F(iV-2,/) 
jj' 

where F denotes the ground band normalization ([T3|) for — 1 and N — 2 bosons. 

The expectation value of a scalar operator O in the 7 band (pOD, with angular 
momentum projection, is given by 

or _(_ 1 /• ( 

{0),,L = _ L) J df3smpdUf3)[{0\b'''%de-^^'^y{bY-'bl\0) 

+2^^{0\b''-%lde-'''''y{bY''bl\0) 
+e^(0|6^-26?6e-^^^-(6t)^-2(6t)2|o)l (33) 



where M{(j)^, L) is the normalization for the 7 band, obtained from eq. ( P3[ ) using 
the identity operator for O. The contribution from the orthogonality terms to the 
band energies are of the order and therefore they were ignored in the original 

papers |]ll| . In description of high-spin states, however, these terms make essential 



contributions and they have to be included in the calculations. Each contraction 
of the intrinsic boson operators in ( ^31) leads to projected single-boson overlaps of 
the form XimXim'd}mm' ■ The resulting Wigner (i-functions are coupled to a final d- 
function to perform the (3 integral. This process leads to rather long expressions 
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for the orthogonality terms. In order to reduce their size, we introduce a compact 
notation for the recouphng coefficients as follows 

R2{jmm', Inn'; I) = XjmXjm'XinXin'{jniln\Im + n){jm'ln'\Im' + n'), 
R^ljmm', Inn', kfj,fi'; I, J) = R2{jmm' , Inn'; I)xk^Xk^i{I'm + nk^\Jm + n + fj,) 

{Im' + n'kjj'lJm' + n' + jj,'), 
R^ljmm', Inn', kjin' , fcVz/'; I, I', J) = R2{jmm' , Inn'; I)R2{kfifi', k'^v'; I') 

{Im + ra/'yU + v\Jm + n + yU + z/) 
{Im' + n7 V + v' I Jm' + n' + ^i' + v) . (34) 

Higher recoupling coefficients [R^,R^ are defined similarly. Using this notation, 
the reduced normalization for the 7 band, F^{N,V) = A/'(0^, L)/(2L + 1), can be 
written as 

F,{N,L) = Y^{L2j-2\I0)' 

X \x%F{N -l,I) + {N-l)J2 ^2(/'20, /02; j)F{N - 2, 1) 



N -I 



2Y,R2{l'2l,mi;j)F{N-2,I) 

- w 

+ {N -2) J2 R3{kW,k'10,m;l',j)F{N -3,1) 

kk'lV 

Y.R2{l'll,lll;3)F{N-2,I) 

w 

+2{N -2) J2 R3{kW,k'01,lll;r,j)F{N -3,1) 



kk'W 



+ -{N -2){N -3) J2 RiiklO,k'lO,m,l'01;k",l",j) 

^ kk'k"ll'l" 

xF{N-4,I) |. (35) 

Eq. (^) expresses the 7-band normalization in terms of the ground-band normal- 
ization ([131) , and it can be evaluated to any order in using Mathematica. 

The expectation value of the number operator in the 7 band can be calculated 
similarly, giving 



(ni) 



1,L 
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+ {N - 1) ^ i?2(/'22, 100; j) + 2/?2(/'20, 102;j)]F{N - 2, /) 
^{N -\){N -2) J2 R3{k20, k' 02, 100; I', 3)F{N -3,1) 



kk'V 



2 ^ i?2(/'21, lOl; j) + i?2(/'01, /21; j)]F{N - 2, /) 

+ (A^ - 2) ^ (2R3{k2l, k'Ol, 100; j) + 2i?3(A;20, A;'01, /Ol; j) 



+i?3(fcl0, k'lO, 102; r, j) j F(A^ - 3, 1) 
+ {N-2){N-3) RA{kl0,k'10,l'02,l00;k",l",j)F{N -4,1) 

kk'k"l'l" 



A^- 1 



2j2R2ii'n,in;j)FiN-2,i) 
I' 

+ {N - 2) ^ (Rsikll, k'll, 100; l',j) + 4i?3(A;ll, k'lO, 101; l',j) 



kk'V 



+2R3{klO, k'Ol, 111; l',j))F{N - 3, /) 



+2(A^-2)(iV-3) (R4{kll,k'lO,l'01,m;k",l",j) 

kk'k"Vl" ^ 



+Ri{klO, k'lO, I'Ol, 101; k" , f,]) ]F{N - 4, /) 



(36) 



It can be easily checked that the condition, I]«(^i)7,L = is satisfied by (|36|) . The 
expectation value of a general two-body interaction in the 7 band is given by 



21 + 1 



(Ar -2)(2A: + 1) 



X (iV-1) 



N-2 



(P2002 + Poo22)F{N - 2, 1)6.. 



iO 



+ (^4Poooo + 4:Xi2XiP2ooo^F{N - 3, 1) 

+ {N -3)Y R2{k20, k'02; t)PooooF{N - 4, /) 



kk' 



A^-2 



P2ioAoF{N-2,I) 



+ (2Xi2XiiPoooi + Xi2XiPoioi + 2XiiXi{P2001 + Po02l)^ F{N - 3, /) 

+ (iV - 3)J2(R2{k21, k'Ol; ^)Poooo + 2i?2(A;20, k'Ol; z)Poooi 

kk' ^ 
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+i?2(A;10, k'lO; z)Pooo2 F{N - 4, /) 



+l(Ar-3)(iV-4) J2 R3{klO,k'lO,t'02;k",i)PooooF{N - 5,1) 

^ kk'k"i' 



-PuiiF{N-2,I) 



+ ^'^ 


■ 1 


iV- 1 


iiV- 2 



2 x2,(Piooi + Pooii) + 2xax,Pioii F(Ar - 3, J) 



+ (iV - 3)Y.i-R2{kll, k'll; z)Poooo + 4^2(^:11, k'lO; z)Poooi 

+2P2(A:10, k'Ol; t){P,ooi + Pooii) + R2ikW, k'lO; t)Poioi )F{N - 4, /) 
+ {N-3){N-A) J2 (R3{kn,k'W,m;k",t)Poooo 



kk'k"i' 



+2P3(A;10, k'lO, i'Ol; k" , i)Poooi j F{N - 5, J) 

+-(A^-3)(A^-4)(iV-5) RA{klO,k'10,t'01,i"01;k",l",t) 
^ kk'k"i'i"i" 



Here we have introduced the compact notation for the two-boson m.e. 



(37) 



{j'm'jm\Jm + m'){lnl'n'\Jn + n') 
{Jm + m'i2-m- m'\I'2) {Jn + n i2 - n - n'\I'2). (38) 

The dummy summation indices I, I', J, i, I' in P are suppressed for convenience. 
The first term in eq. ( p7| ) is the effective one-body term that arises from normal 
ordering of the boson operators in the multipole interaction, and it is expressed 
using eq. (P^l). 

Eqs. ( ^|37D are the counterparts of eqs. ([TTl , p!7D for the ground band and can 
be fed directly into Mathematica for evaluation. As one can surmise from a cursory 
comparison of the parent equations, the resulting third layer expressions are pages 
long. They are not as accurate as the ground band results, presumably requiring 
inclusion of even higher order terms. For these reasons, we have opted for a numerical 
evaluation of eqs. (3^; 37) in the applications. Such a calculation includes all orders 
in and hence provides more reliable results for the 7-band energies. 
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4.4 E2 transitions 



Description of the yrast E2 transitions at high-spins is one of the main aims of this 
work. Therefore, we present a brief review of the currently available results 
for E2 m.e. and discuss their extensions to higher orders. A comprehensive study 
of the E2 transitions among the ground, 7 and (3 bands was given previously |j23|. 



The first layer m.e. obtained in for the yrast E2 transitions appears to work 
rather well even at high-spins Inclusion of the d-boson energy leads to some 
deterioration at high-spins, which can be rectified by incorporating the higher order 
terms in the expansion. 

The ground band m.e. of the quadrupole operator is given by 

{L' II Q \\ L) = L[AF{N, L')F{N, L)]"^/^ ^(^^2 - M|L'0) 

■ci/3sin/54,o(0|6^Q_A/e-'^^''(6t)^|0), (39) 

where L = [2L + 1]^/^. As before, this can be reduced to the form 

{L' II Q\\L)= [^(^^^^(^j^^^)]i/2E?.'^.^^ 

x{jOL'0\JO){mO\JO)\ I f ■!,]f{N-1,J), (40) 



which can be evaluated to any order using Mathematica. Because of the tensor 
nature of the E2 operator, however, the resulting third layer expressions are much 
more complicated than those for the Hamiltonian. In contrast, because there is 
no variation involved, numerical evaluation of eq. (|^) is straightforward and is 
preferred over the lengthy algebraic forms in the following. 

4.5 Comparison with the exact results 

Before applying the 1/N expansion results, we compare them with those obtained 
from an exact diagonalization of the Hamiltonian ||l^. Of necessity, the boson 



number is fixed at = 10. The Hamiltonian parameters are as in sect. 3 {k,2 = —20 
keV, K4 = 0, q = 0.5, rjd = 1.5, rjg = 4.5), except where noted. Fig. 3 shows a 
comparison of the ground band energies normalized with L = L{L + 1) so that 
they all have the same energy scale. Fig. 3a studies the convergence of the 1/A^ 
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results obtained with the YAP procedure. The second layer results (dotted 
line) rapidly diverge from the exact energies (circles) for spins L > 2N, and hence 
are not reliable in applications to high-spin states. The third layer results (dashed 
line), on the other hand, track the exact energies within a few percent up to the 
maximum spin L = 4N. Using the numerical technique described in Appendix 
A, one can evaluate the expansion to all orders (solid line) which exhibits an 
almost perfect agreement with the exact energies. This study demonstrates that 
the third layer expansion results are both necessary and sufficient for a reliable 
description of high-spin sates. Fig. 3b shows the effect of the one-body energies on 
the accuracy of the third layer results. The top line is as in fig. 3a, the middle one 
compares the exact and results for rj^ = 0, r]g = 4.5, and the bottom one for 
^7^ = 0, ?7g = 0. It is seen that the agreement for a pure quadrupole Hamiltonian 
is excellent at all spins, while the addition of g-boson energy leads to a few percent 
deviation at very high-spins. In a typical situation with d-boson energy, this few 
percent deviation starts occurring at medium high-spins. 

In fig. 4a, we present a similar study for the gamma band energies. The av- 
erage behaviour is well reproduced by the expansion results but staggering is 
underestimated. This happens because staggering is caused mainly by band mixing 
between the ground and 7 bands which is not included in the present calculations 
(note that the odd-spin levels, which are not affected by band mixing, are very well 
reproduced). We have not attempted to include band mixing effects here because 



they are strongly suppressed for the larger values used in deformed nuclei p3| , 
and hence they can be ignored for the purposes of this work. Finally, in fig. 4b, 
we compare the 1/A^ results for the yrast E2 transition m.e. with the exact ones 
(circles). The dashed line shows the first layer result obtained from eq. (^2]) which is 
accurate to a few percent for L < 2N, but progressively gets worse with increasing 
spin. Solid line shows the numerical evaluation of eq. (0), which is complete to all 
orders in 1/A^. The agreement with the exact results becomes almost perfect in this 
case, including the highest spins which are dominated by g bosons. 
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5 Systematic studies 



Application of the IBM to deformed nuclei has been criticized on a number of fronts 
ranging from the energy scale of 7 excitation modes to the spin dependence of 
MOT and E2 transitions 0]. Most of these problems have been attributed to the 
truncation of the model space to s and d bosons but due to lack of adequate tools they 
have not been properly addressed in the framework of the sdg-IBM. The analytic 
formulas obtained in the last section have the advantage that one can easily perform 
systematic studies of key physical quantities and obtain useful insights on the effect 
of various parameters. In this section, we present such a study that will shed light 
on the above problems and suggest more appropriate Hamiltonians for description 
of deformed nuclei. As mentioned in sect. the minimal sdg-IBM Hamiltonian 
with quadrupole interaction and g-boson energy has previously been applied to 
deformed nuclei and found to be inadequate on most counts |]12|, |T^. Inclusion of 
the d-boson energy and hexadecapole interaction appear to be necessary to improve 
this situation. Here we focus on the effect of these extra terms on the minimal 
Hamiltonian, especially with regard to high-spin sates. 

To simplify the discussion, we rewrite the ground and 7 band energies as 

where the coefficients A„ can be read from the respective energy expressions. Eq. ( PID 
is the familiar rotational expansion of the level energies used in the geometrical model 



2^ . The difference between the two models is that in the IBM the coefficients A„ 
follow from an underlying Hamiltonian (which is used in describing other properties) 
whereas in the geometrical model they are directly extracted from the data. The 
MOI problem raised in Ref. |^ refers to the fact that i) the Ai coefficient gets a 
substantial contribution from the dipole interaction, L ■ L, which has no dynamical 
content, ii) the A2 coefficient is much smaller than the experimental values, and iii) 
the variation in Ai among different bands can not be described. All three problems 
are in fact interrelated. Although the second can be resolved by renormalizing the 
moment of inertia at high-spins (e.g. by modifying L ■ L L ■ L/{1 + fL ■ L) 



|25[|), such modifications are purely kinematical in origin and do not address the 
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dynamical problem. Further quantities of interest in the study of high-spin states 
are the yrast E2 transitions. For systematics, it is sufficient to consider the first 
layer expansion result which has the generic form 

(L - 2 II T{E2) II L) = e2NL{L0 20|L - 2 0) [rrii + m2L{L - 1)] (42) 

where the coefficients m„ are given in ref. The first term in (E2p gives the 



familiar rigid-rotor result. The second term is negative and is responsible for the 
falloffs predicted in E2 transitions. 

In presenting systematics, we find it convenient to use ratios which eliminate the 
undesired effects of the scale parameters k,2 and N. The energy scale can be fixed, 
for example, by fitting K2 to the excitation energy of the 7 band, E^. We discuss 
five such ratios as a function of q for various values of (a) rjd = ed/NK2 and (b) 
C4 = K4/k2- The parameter q is varied from 0-1 which covers the whole range of 
the quadrupole operator from the 7-unstable to the SU(3) limit, rid is varied from 
0-2 in 10 equal steps, and ^4 from 0-0.5 in 5 equal steps, which cover the range of 
values used in the applications. (Negative values of (4 on the whole are found to 
have an adverse effect and hence are not considered.) In the rjd systematics study, 
C4 = is used as its precise value does not have much influence on the results. In 
the ^4 study, however, the choice of rj^ does have an impact, and we adapt 77^ = 1.5 
which is the average value used in the applications. The g-boson energy has been 



found to have negligible effect |]14|, and therefore is not varied here but flxed at 
r]g = 4.5. Below, we comment on the behaviour of each ratio and contrast them 
with the experimental data. For reference, we note that q assumes values around 
~ 0.5 in the rare-earth nuclei and ~ 0.7 in the actinides. 

1) E^/NXgi (flg. 5): This ratio relates the energy scales of the 7 and ground 
bands, and its mismatch with experiment has been a source of criticism ||^. It is 
around 4-5 in the rare-earth region and increases to 8-9 in the actinides. From flg. 5a 
it is seen that the CQF with rjd = overestimates it by about a factor of 2 in both 
regions. It decreases rapidly with rja though, and through a judicious use of the 
d-boson energy, it should be possible to describe this ratio (and hence the MOI) 
without appealing to the L ■ L term. Fig. 5b shows a similar study on the effect of 
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the hexadecapole interaction which is seen to be going in the right direction but is 
too small to have any impact. 

2) N'^Xg2/Xgi (fig. 6): This ratio measures the deviation from the rigid rotor 
behaviour due to loss of pairing. It ranges from about -0.2 in the rare-earth region 
to -0.1 in the actinides. The CQF with rjd = gives values an order of magnitude 
smaller (fig. 6a) and hence signally fails in accounting for the spin dependence of 
MOI as first pointed out in ref. However, this ratio is very sensitive to the rja 
values and the experimental range can be easily attained by including the d-boson 
energy in the Hamiltonian. From fig. 6b, the hexadecapole interaction is seen to have 
a coherent effect in further reducing this ratio away from the rigid rotor behaviour. 

3) N'^\^2/\i (fig- 7): An identical study for the 7 band indicates broadly sim- 
ilar but somewhat larger effects of the d-boson energy on the behaviour of 7 band 
MOI (fig. 7a). A softer MOI in the 7 band is in line with data in most deformed 
nuclei though there are a few exceptions as will be seen in the applications. The 
hexadecapole interaction has an opposite effect (fig. 7b) which reduces the difference 
between the ground and 7 band MOI caused by the d-boson energy. 

4) X-yi/Xgi (fig. 8): This ratio compares the MOI of ground and 7 bands which 
fiuctuates within a band of ±10% across the deformed nuclei. The earlier IBM 
calculations gave results near one and could not accommodate such fiuctuations. 
Fig. 8a-b shows that inclusion of the d-boson energy can increase this ratio by up 
to 20-30%, while the hexadecapole interaction can reduce it significantly (up to 
20-30%) thereby covering the whole range of fiuctuations. 

5) N'^m2/mi (fig. 9): As there is no boson cutoff effect, experimentally this ratio 
is consistent with zero. The CQF with rjd = leads to values similar to the sd- 
IBM (fig. 9a), hence despite the addition of g bosons, it would suffer from the same 
boson cutoff problem. Introduction of the d-boson energy, however, reduces it sub- 
stantially, becoming more in line with experiments. The effect of the hexadecapole 
interaction on this ratio (fig. 9b) is similar to fig. 5b; it is positive but comparatively 
too small to make a difference. 

Another ratio, namely, N^Xg^/Xgi is also of interest especially at very high-spins 
(L = 20 — 30) where the cubic term in eq. ( ^T]) plays an important role |T^. It 



23 



exhibits a similar dependence on 77^ and ^4 as N'^\g2/^gi in fig. 6, so it is not discussed 
further here. 

The d-boson energy has been mostly neglected in studies of deformed nuclei, 
presumably due to the success of the CQF with Ed = in explaining the energy and 
E2 transition systematics of low-lying states In fact, for small values {rid ~ 1), 
its effect on low-lying states is negligible and it is not really needed in their descrip- 
tion |]2B[. The CQF, however, basically leads to a rigid MOI and can not explain 
either its spin dependence or its variation among different bands. The obvious way 
towards a softer energy surface is to include the d-boson energy in the Hamiltonian 
which is seen to vastly improve the description of the spin-dependent terms in the 
level energies and E2 transitions. The hexadecapole interaction performs a similar 
function but has a much smaller effect. The exception is, of course, variations in 
MOI which could not be reproduced without the hexadecapole interaction. 

6 Applications to deformed nuclei 

In the light of the systematic trends discussed above, we carry out fits to the rare- 
earth nuclei i58-i62j3y^ 164-168^^, les-iTGyb, iTO-irsjjf and the actinides 228-232^^^^ 
234-238-[j^ The isotopes chosen are all well deformed rotors with energy ratio E4/E2 
close to 3.3. We have excluded those exhibiting backbending as their proper de- 
scription requires inclusion of two-quasiparticle states in the model space. While we 
mainly focus on the description of high-spin states, which has not been done before, 
we also consider a selected set of low-lying bands. This is important in properly con- 
straining the model parameters so that the results obtained are valid in a broader 
sense and not just for a small subset of observables. The sdg-IBM parameters used 
in the fits are listed in table 1. Each parameter is particularly sensitive to a certain 
set of observables which simplifies the fitting process. For example, q is determined 
from interband E2 transitions, rjd from the spin dependence of MOI and E^^ (cf. figs. 
5-6), rig from £'3+, E4+, (4 from MOI variation (cf. fig. 8), and finally K2 from the 
overall energy scale of the spectrum. The parameters are either constant in a given 
isotope chain or change smoothly in accordance with the vibration-rotation shape 
transition, e.g., rjd decreases with increasing as the nuclei considered become more 
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rotational. 

The representative observables chosen to describe the low-lying band structure 
are the band excitation energies Ep, E^, E3+, E4+ (table 2), the interband E2 m.e. 
for 2/3 ,y Og transitions (table 3), and the E4 m.e. for 4^,3+ ,4+ Og transitions 
(table 4). Note that the E4 m.e. are normalized with the ground transition, so that 
an effective E4 charge is not needed in table 1. With a few exceptions to be discussed 
below, the general trends of the (3 and 7 band systematics are well reproduced by 
the calculations. The sudden fluctuations seen in some of the band-head energies 
(table 2) can be accommodated by a careful tuning of the parameters. Since our 
aim here is to delineate the systematic features of high-spin states, rather than to 
obtain refined fits to individual nuclei, we have not attempted such an improvement. 
Description of interband E2 transitions is one of the strong points of the IBM, and 
as can be seen from table 3, they are very well reproduced using almost constant 
q values. The parameters in the E4 operator are determined from the conditions 
(^, hence the E4 m.e. ratios presented in table 4 are parameter free predictions 
of the model. Again the overall agreement with the data is reasonable which gives 
confidence on the choice of the E4 operator. The quality of agreement obtained 
in tables 2-4 indicates that the limited set of sdg-IBM parameters (table 1) can 
describe the basic features of the low-lying bands in deformed nuclei. 

In the study of high-spin sates, we include the level energies for the ground and 
7 bands, and the yrast E2 transitions for each set of isotopes (figs. 10-27). We 
first comment on their general features. In all cases, the MOI strongly deviates 
from the rigid rotor behaviour which would be represented by a horizontal line in 
the figures. Further, this deviation is not linear but curves up with increasing spin 
underscoring the importance of the cubic term in eq. (|4lD . Note that because of the 
ample data available, these features are most clear in the ground bands and to a 
lesser extent in the 7 bands. The yrast E2 m.e., on the other hand, follow closely 
the rigid rotor values with no sign of a boson cutoff effect. As emphasized in sect. 
5, these properties can be explained in the IBM by including the d-boson energy in 
the Hamiltonian. Below we comment on the specific features of each isotope chain. 

1) i58-i62£)-y (ggs, 10-12): Among the deformed nuclei considered in this work. 
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^^^Dy, together with ^ Hf, exhibit the largest changes in MOI. These nuclei have 
the lowest boson numbers among the rare-earth set and are clearly influenced by 
the vibration-rotation phase transition as indicated by the larger rjd values used. 
The ground band energies (flg. 10) are well described with relative errors of about 
1-2%. The trend in 7 band energies (flg. 11) is similarly reproduced (note the 
different scales in flgs. 10 and 11). The slight overprediction of energies here can 
be improved by flne tuning the hexadecapole interaction (cf. flg. 8b). The yrast E2 
m.e. have been a sore point in applications of the sd-IBM to high-spin states due 
to boson cutoff. For example, in ^^'^Dy, the sd-IBM would predict band termination 
at L = 26 which is not seen in the data (flg. 12a). This problem has been resolved 
in the present sdg-IBM calculations which account for the yrast E2 data very well 
(flg. 12). A side remark for ^^^Dy is that the band excitation energies in this nucleus 
do not follow the trend of i^s-ieoj^y (table 2), hence it requires individual attention 
for a better description. 

2) 164-168 j^j^ (flgs. 13-15): The Er isotopes, and in particular ^^^Er, are the 
exceptional cases mentioned above for which a consistent description of the data 
could not be obtained with our limited set of parameters. While the spin dependence 
of the ground and 7 band MOI (flgs. 13-14) and the E2 m.e. (flg. 15, table 3) are 
well described, the band excitation energies are overpredicted (table 2) and the E4 
m.e. are rather poor (table 4). The problem stems from the fact that, among all 
the deformed nuclei considered in this study, the Er isotopes have the lowest lying 
7 bands and the most rigid MOI. As seen from flgs. 6-7, these two quantities are 
correlated in the present parametrization, so that a lower 7 band obliges a softer 
MOI (cf. flgs. 6-7). Thus a proper description of the Er isotopes requires extension 
of the Hamiltonian (^ and/or relaxation of the constraints on the quadrupole and 
hexadecapole parameters. For example, in a detailed study of ^^^Ei in the sdg- 
IBM, 14 parameters were employed 0. Here we will be content with exposing the 
exceptional nature of the Er isotopes and leave their detailed investigation for future 
work. 

3) 168-176 Y]-) (flgs. 16-18): The Yb isotopes are uniformly well described and 
require little comment. In contrast to Er, the 7 band energies in the Yb isotopes are 
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higher which are well correlated with their relatively stiff MOI (fig. 16). One point 
worthwhile to make is that the 7 band MOI is higher than that of the ground band 
which could not be explained without the hexadecapole interaction. 

4) I'i'o-iTSjjf (figs. 19-21): In the Hf isotopes, the 7 band comes down but the 
MOI are softer, and hence the correlation between the two quantities is preserved 
(fig. 19). The staggering observed in the 7 bands (fig. 20) requires inclusion of band 
mixing effects for a better description. Otherwise the data are well reproduced by 
the calculations. 

5) 228-232rj]-^ (figs. 22-24): Although boson numbers are lower in the actinide 
nuclei considered here, they exhibit characteristics of well deformed nuclei. The 
MOI in actinides are typically twice as large as those in rare-earths requiring smaller 
K,2 values. The high-spin data are scarce in 228-230'j]-^ ^^^Th, where data up 
to spin L = 30 are available, an excellent description is obtained. One interesting 
feature of the 7 band MOI in ^^^Th is that it is larger and stiffer compared to the 
ground band. Both of these features require a large dose of hexadecapole interaction 
for explanation as remarked in the systematic studies. 

6) 234-238U (^ggg_ 25-27): The most extensive high-spin data are available for the 
yrast bands in the U isotopes which are well described by the present calculations. 
The yrast E2 m.e. in the U isotopes (and ^^^Th) were measured to check the boson 
cutoff predictions of the sd-IBM, i.e. E2 m.e. vanish at L = 2N. As seen in fig. 27, 
the E2 data show no sign of falloff which provides one of the strongest motivations 
for inclusion of g bosons. At the highest spins, the sdg-IBM calculations appear to 
underpredict the E2 measurements. We emphasize that this is not due to any boson 
cutoff effect but rather due to deviation of the data from the rigid rotor values. That 
is, most models would have difficulty in explaining these E2 transition m.e. which 
are larger than the rigid rotor values (see, for example, p3[ ). 

7 Summary and conclusions 

In this paper, we have presented a systematic description of high-spin states in de- 
formed nuclei within the framework of the sdg-IBM. Such a description has been 
long overdue but has not been performed earlier due to the technical difficulties in 
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diagonalizing the sdg-IBM Hamiltonians. We have shown that the 1/N expansion 
formahsm, extended to higher orders using computer algebra, provides a viable al- 
ternative for this purpose. Systematic studies of the model parameters using the 
expansion formulas have indicated that some of the long standing problems as- 
sociated with the description of MOI and E2 transitions in the IBM can be resolved 
by including the d-boson energy in the Hamiltonian. The hexadecapole interaction 
has a minor effect on the ground band but could play a decisive role on 7 and 
other excited band properties. This feature of the hexadecapole interaction does 
not appear to be well appreciated in literature. For example, the recently observed 
staggering effect in some superdeformed bands have been attributed to the hexade- 
capole degrees of freedom. If this is true, then it would have profound effects on the 
neighbouring non-yrast bands. 

The apphcation of the sdg-IBM Hamiltonian consisting of single-boson ener- 
gies, and quadrupole and hexadecapole interactions (with constrained parameters) 
resulted in a mostly uniform and successful description of both low- lying band struc- 
tures and high-spin states across the rare-earth and actinide regions. In the past, 
many experimental results on high-spin states were compared with the sd-lBM cal- 
culations with negative connotations. This was presumably due to the lack of the 
sdg-IBM calculations. We hope that the extensive sdg-IBM results presented in this 
work will help remedying this situation. 

Our limited parametrization did not work as well in the case of the Er isotopes 
which appear to have rather exceptional properties as exposed in this systematic 
study. We are planning to carry out a more detailed study of the Er isotopes without 
imposing the constraints used in this paper. It is amusing to note that ^^^Er, which 
has been used as a benchmark case in tests of phenomenological models, may turn 
out to be the most exceptional of all deformed nuclei. 

Appendix A 

In cases where the higher order 1/N expressions are too complicated to be useful, 
numerical evaluation of the 1/N expansion formulas may provide a more practical 
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alternative. The reduced normalization integral ([T2|) can be written explicitly as 



F{N,L)= [ 
Jo 



1 



(43) 



where Pl{z) is a Legendre polynomial. In the sd-lBM, the term in brackets is given 
by (co + C2Z^)'^ , and in the sdg-IBM by (cq + C2Z^ + c^z'^)^ where q are constants 
involving the mean fields. By expanding the binomial or trinomial, and using 



hence {H) are calculated exactly to all orders in The variational problem 

is solved numerically using the exact m.e. (|ll],|l3) and the simplex method. The 
numerical evaluation can be used, for example, to check the convergence of 
expressions and determine whether a calculation to a particular order is accurate 
enough. 




(44) 



eq. ( P3| ) can be integrated term by term. This has the advantage that F{N, L) and 
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Figure captions 

Fig. 1. Effect of the basis space truncation on a) band excitation energies b) E2 
transitions, and c) E4 transitions. The maximum number of g bosons, n^max, al- 
lowed in the basis space is increased from 1 to the maximum of = 10. Parameters 
of the sdg-IBM Hamiltonian are given in the text. 

Fig. 2. Effect of the basis space truncation on the ground band a) excitation ener- 
gies, and b) E2 transitions. 

Fig. 3. Comparison of the ground band energies obtained from the expansion 
with the exact diagonalization results (circles). In (a) different lines refer to the sec- 
ond layer calculation (dotted line), the third layer (dashed line), and the numerical 
one to all orders (solid line). In (b) the lines correspond to the third layer results 
obtained with r/^ = 1.5, rjg = 4.5 (top), rj^ = 0, r]g = 4.5 (middle), rj^ = 0, rjg = 
(bottom). 

Fig. 4. (a) Comparison of the 7 band energies obtained from the 1/A^ expansion 
with the exact diagonalization results (circles), (b) Comparison of the yrast E2 
transition m.e. obtained from eq. (|4^ ) (dashed line) with the exact diagonalization 
results (circles). The solid line shows the numerical evaluation of the m.e. to all 
orders. 

Fig. 5. The effect of (a) the d-boson energy and (b) the hexadecapole interaction 
on the ratio E^/NXgi which relates the energy scales of the 7 and ground bands. 

Fig. 6. Same as fig. 5 but for the ratio N'^Xg2/Xgi which measures the deviation of 
the ground band moment of inertia from the rigid rotor behaviour. 

Fig. 7. Same as fig. 6 but for the 7 band. 

Fig. 8. Same as fig. 5 but for the ratio A^i/Agi which compares the moment of 
inertia of the ground and 7 bands. 

Fig. 9. Same as fig. 5 but for the ratio m2/mi which measures the boson cutoff 
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effect. 



Fig. 10. Comparison of the experimental (circles) and calculated (solid lines) ground 
band energies EgL/L{L + 1) (in keV) in ise-ieoj^y^ rj.^^ ^g^^g^ ^^^^ l^, |28|, |29 . 



Fig. 11. Same as fig. 10 but for the 7 bands. 



Fig. 12. Comparison of the experimental (circles) and calculated (solid lines) yrast 
E2 transitions in ise-ieo^y rj.^^ ^g^^g^ ^^.^ ^^.^^ 0, ^ ^ . 



Fig. 13. Same as fig. 10 but in Er. The data are from |30, 31, p2| 



Fig. 14. Same as fig. 11 but in ° ^°^Er. The data are from |30, 31, p2^ 



Fig. 15. Same as fig. 12 but in ^^^Er. The data are from ||30|, |31|, B2 



Fig. 16. Same as fig. 10 but in ^^^-^^^Yb. The data are from ^ |3|, ^ |36 



Fig. 17. Same as fig. 11 but in ^^^-^^^Yb. The data are from [32, 33, 34, 35, 36 



Fig. 18. Same as fig. 12 but in les-ireyb. The data are from H, H, |4H|, |3|] . 



Fig. 19. Same as fig. 10 but in ^^'^'^^Ri. The data are from ||, ||, ||, 



Fig. 20. Same as fig. 11 but in Hf. The data are from |RB|, B3, 



Fig. 21. Same as fig. 12 but in 



i70-i78Hf. The data are from |3|, p|, p5l ||, |7 



Fig. 22. Same as fig. 10 but in ^^«-^^^Th. The data are from ||, ||, |g . 



Fig. 23. Same as fig. 11 but in 228-232rjj^_ rj.^^ ^g^^g^ g^^^ ^^^^ ^ 



Fig. 24. Same as fig. 12 but in 228-232r,nj^_ ^^le data are from m, M, p 



Fig. 25. Same as fig. 10 but in 234-238u. The data are from ET], EOl |2 . 



Fig. 26. Same as fig. 11 but in 234-238u. The data are from gl], |0|, |2 



Fig. 27. Same as fig. 12 but in 234-238u. The data are from gl], H, |2| 
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Table 1: Parameters used in the sdg-lBM calculations. K2 is in keV and 62 in eb. 



AT 1 

iNucleus 


A T 

N 


K2 


C4 


5 


rid 




62 


158 


13 


19.8 


0.30 


0.50 


1.90 


5.0 


0.13 


160 


14 


19.9 


0.30 


0.50 


1.77 


4.6 


0.13 




15 


19.3 


0.35 


0.50 


1.60 


4.5 


0.13 


164 


14 


22.0 


0.40 


0.50 


1.50 


5.0 


0.14 


166 


15 


21.2 


0.40 


0.50 


1.42 


4.7 


0.13 


168 


16 


23.1 


0.40 


0.50 


1.22 


4.6 


0.13 


168Yb 


14 


19.7 


0.35 


0.50 


1.68 


5.0 


0.14 


170Yb 


15 


20.6 


0.35 


0.50 


1.63 


4.7 


0.13 


i^^Yb 


16 


20.5 


0.25 


0.60 


1.78 


5.0 


0.13 


174Yb 


17 


22.2 


0.25 


0.60 


1.66 


4.1 


0.12 


176Yb 


16 


20.1 


0.35 


0.60 


1.83 


5.3 


0.12 


170Hf 


13 


19.1 


0.10 


0.50 


2.04 


4.9 


0.14 


172Hf 


14 


19.5 


0.10 


0.50 


1.99 


4.5 


0.13 


174Hf 


15 


20.7 


0.10 


0.50 


1.91 


4.1 


0.13 


176Hf 


16 


22.0 


0.10 


0.50 


1.80 


3.5 


0.13 


178Hf 


15 


21.8 


0.10 


0.50 


1.83 


4.1 


0.13 


228xh 


10 


19.7 


0.30 


0.68 


1.60 


3.0 


0.20 


230 rpj^ 


11 


15.5 


0.40 


0.68 


1.59 


4.4 


0.20 


232 rpj^ 


12 


14.6 


0.40 


0.68 


1.58 


4.4 


0.20 


234U 


13 


14.9 


0.20 


0.70 


1.62 


3.5 


0.18 


236 u 


14 


16.4 


0.20 


0.70 


1.57 


3.2 


0.17 


238^: 


15 


17.7 


0.20 


0.70 


1.56 


2.8 


0.17 
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Table 2: Comparison of the /5, 7, and K = S"*", 4"*" (single-phonon) band energies 
(in keV) with the sdg-IBM calculations in the rare-earth and actinide regions. The 

data are from Nucl. Data Sheets. 

Eg E^+ E^+ 



Nucleus 


Cal. 


Exp. 


Cal. 


Exp. 


Cal. 


Exp. 


Cal. 


Exp. 


158j3y 


916 


991 


965 


946 


1461 




1935 


1895 


leo^y 


1204 


1275 


998 


966 


1512 




2085 






1284 


1205 


1048 


888 


1617 




2201 


1536 


J_J1 




1246 


1 1 45 


860 


1 634 


1702 


2003 






1275 


1460 


1138 


786 


1704 




2376 




168 


1590 


1217 


1345 


821 


1892 


1654 


2483 


2238 


168Yb 


1085 


1156 


1028 


984 


1470 


1452 


1706 




170Yb 


1218 


1069 


1138 


1145 


1501 




1857 


1408 


i72Yb 


1345 


1043 


1385 


1466 


1799 


1663 


2403 


2073 


i^^Yb 


1589 


1487 


1554 


1634 


1810 




2545 




176Yb 


1329 


1779 


1324 


1261 


1685 




2407 




170Hf 


868 


880 


985 


961 


1482 




1528 




172Hf 


982 


871 


1016 


1075 


1530 




1604 




174Hf 


1108 


827 


1137 


1227 


1547 


1303 


1669 




176Hf 


1298 


1150 


1278 


1341 


1671 


1578 


1840 




178Hf 


1231 


1199 


1202 


1175 


1710 


1728 


1809 


1848 


228xh 


791 


832 


846 


969 


1082 




1458 




230 


734 


635 


792 


781 


1004 




1519 




232Th 


769 


730 


813 


785 


1041 




1574 




234u 


825 


809 


871 


927 


1253 


1496 


1566 


1723 


236^ 


986 


919 


1002 


958 


1325 




1799 




238^ 


1080 


993 


1086 


1060 


1328 


1059 


1965 
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Table 3: Comparison of the interband E2 transitions (in eb) with the sdg-IBM 
calculations in the rare-earth and actinide regions. The data are from Nucl. Data 

Sheets. 

(2^||T(i?2)||0,) (2,||T(E2)||0,) 



Nucleus 


Cal. 


Exp. 


Cal. 


Exp. 


158 


0.21 


0.23 ± 0.02 


0.41 


0.39 ± 0.04 


leo^y 


0.19 




0.34 


0.27 ± 0.04 




0.17 




0.35 


0.35 ± 0.02 




17 




34 


37 + 02 




0.16 




0.41 


0.39 ± 0.02 




0.15 


<0.03 


0.34 


0.36 ± 0.01 


168Yb 


0.20 


0.22 ± 0.01 


0.35 


0.36 ± 0.04 


170Yb 


0.19 


0.17 ± 0.02 


0.32 


0.28 ± 0.03 


i^^Yb 


0.16 


0.09 ± 0.01 


0.26 


0.21 ± 0.03 


i^^Yb 


0.14 




0.24 


0.22 ± 0.03 


176Yb 


0.14 




0.27 


0.23 ± 0.03 


170Hf 


0.23 




0.36 




172Hf 


0.19 




0.31 




174Hf 


0.20 




0.33 


0.37 ± 0.04 


176Hf 


0.18 


0.21 ± 0.02 


0.34 


0.35 ± 0.01 


178Hf 


0.21 




0.32 


0.34 ± 0.02 


228 


0.18 




0.28 




230 


0.20 


0.21 ± 0.05 


0.34 


0.35 ± 0.06 


232Th 


0.23 


0.31 ± 0.07 


0.38 


0.36 ± 0.04 


234U 


0.21 


<0.24 


0.32 


0.35 ± 0.04 


236 u 


0.20 




0.30 




238^ 


0.19 


0.23 ± 0.03 


0.31 


0.36 ± 0.04 
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Table 4: Comparison of the interband EA transitions, normalized to inband ones, 
with the sdg-IBM calculations in the rare-earth and actinide regions. The data are 
from Nucl. Data Sheets. 

<4^||T(£;4)||03> <43+||T(A'4)||0g-> <4^+ ||T(ji.'4)||03-> 

<4g||T(E4)||0g> <4g\\T{E4)\\0g> <4a||T(E4)||0g> 



Nucleus 


Cal. 


Exp. 


Cal. 


Exp. 


Cal. 


Exp. 


158 J3 


0.49 


56 ± 44 

\J,tJ\J —1— W • -L -L 


1.14 




0.32 




160 J3 


0.50 


63 ± 24 


1.08 




0.27 






0.37 




1.06 




0.28 




164 p 

11/1 


n 44 

U .tit: 




1 .yjo 








166 


0.53 




1.00 




0.18 




168 gj. 


0.46 


1.32 ± 0.72 


0.95 


0.60 ± 0.33 


0.15 


0.25 ± 0.15 


168Yb 


0.50 




0.95 




0.20 




170Yb 


0.60 




0.89 




0.14 




i72Yb 


0.53 


0.20 ± 0.19 


0.82 


0.69 ± 0.57 


0.09 




172Yb 


0.36 


0.20 ± 0.19 


0.78 


0.69 ± 0.57 


0.16 




174Yb 


0.44 




0.73 


0.54 ± 0.29 


0.11 




176Yb 


0.37 




0.81 




0.17 




170Hf 


0.44 




1.17 




0.43 




172Hf 


0.45 




1.12 




0.38 




174Hf 


0.50 




1.04 




0.28 




176Hf 


0.44 




1.01 




0.28 




178Hf 


0.43 




1.08 




0.34 




228 


0.55 




0.62 




0.02 




230 rpj^ 


0.44 




0.71 




0.07 




232 rpj^ 


0.48 




0.68 




0.05 




234U 


0.48 




0.61 




0.02 




238 u 


0.50 




0.58 




0.02 
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